Canted-spin-caused electric dipoles; a symmetry study 
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A pair of magnetic atoms with canted spins S, can give rise to an electric dipole moment A. 
This canting-dependent moment can be written in general as A = Ar x (S x S')-\-BS x S^, where r is 
the relative displacement of the atomic positions. A novel and rigorous general symmetry argument 
shows that A 7^ is allowed no matter how high the symmetry of the atoms plus environment, and 
gives symmetry requirements for B ^ 0. 
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Great recent interest in multiferroic materials, e.g. [l|, 

a i a a a a a a eo, [o, e m b e m 

where magnetic ordering of various sorts induces ferro- 
or ferri-electricity, forces one to understand the mech- 
anism for this surprising, and possibly useful effect. 
One mechanism due to Katsura, Nagaosa and Balatzky 
(KNB) a fllj is derived fundamentally by considering a 
pair of magnetic ions whose spins Si, S2 are constrained 
to be in arbitrary directions. Such a constraint is imag- 
ined to result from exchange fields originating from cur- 
rents or spins outside of the atom-pair considered. E.g, 
the magnetic state in a crystal might be a spiral and the 
atom pair considered would be any neighboring pair par- 
ticipating in the spiral. In ^ El found that the 
electronic density becomes distorted by a combination of 
spin-orbit coupling and interionic electron hopping due 
to the usual kinetic and Coulomb effects, leading to an 
electric dipole moment which is given by 



p = cr x (Si x S2), 



(1) 



where r is the displacement of one atom relative to the 
other, and c is a numerical coefficient. This form was 
derived by others, in the context of an extended system 
with spiral-like magnetic ordering: Derivations based on 
spin-lattice interactions are due to Sergienko and Dagotto 
(SD) [8] and Harris et al [9]. There are also phenomeno- 
logical derivations using symmetry arguments via Lan- 
dau theory Q, and Landau-Ginzberg theory [[3]. A fur- 
ther work relevant here [10] is closely related to the KNB 
mechanism, again involving a pair of atoms; an essential 
difference lies in the fact that the hopping matrix ele- 
ment is that of the spin-orbit coupling (interatomic SO 
coupling), whereas the SO coupling in the KNB approach 
appears only in the matrix involving Wannier functions 
on the same atom (intra-atomic SO coupling). In [10] a 
the relation ([1]) was found, where the assumption was 
made that the spatial symmetry of the situation was the 
symmetry of a pair of points in space, an assumption also 
made in a[ll|; however, in [lO]b, a lower symmetry was 
studied, which led to the possibility of another compo- 
nent of p, namely in the direction Si x S2, thus raising 
the question of the generality of ([1]). This was also raised, 
considering extended systems, in [16], and actually an- 



swered in [13] , where it was shown by an experimental ex- 
ample and a Landau theory that this Si x S2-component 
can exist. Below, the symmetry theory developed here is 
applied to this example, and is contrasted with Landau 
theory in general. 

Because of the complexity of the arguments raised, 
as just described (see also the Reply to [16] by 
Mostovoy [17]), and the intense interest in the general 
topic of multiferroics, we are motivated to give a rigor- 
ous, general, but simple analysis of the conditions under 
which these two types of canted-spin-induced electric po- 
larization can exist. This, which is a pure symmetry anal- 
ysis, is as follows. It is local, in the sense of [5|, llO|, llll |: 
it treats a single pair of ions or atoms. It is closely anal- 
ogous to an argument leading to the stringent conditions 
on the DM vector D (Moriya's rules) imposed by the 
symmetries of the ion-pair plus its surroundings. [18|] We 



FIG. 1: Dots are identical spherical 'atoms' at A and B; x- 
and y-axes point to the right and up respectively. 

begin by reviewing this argument that leads to Moriya's 
rules. [21] One considers the possible existence of a term 
in the energy of the form Edm = D • (Si x S2), where D 
is assumed to be a property of the structure, atom-pair 
plus surroundings and independent of spins (see [18] for 
discussion). Si and S2 are the spins at sites A and B re- 
spectively. One explores the conditions imposed on D by 
possible symmetries of the structure (without spins), i.e. 
rotations which return the two sites plus surroundings 
to itself, with the requirement that Edm be unchanged 
(as a term in a Hamiltonian, it's a scalar under such op- 
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erations). Important is the fact that D is fixed in the 
structure (as seen in Moriya's mathematical expression 
for it), so that D is the same before and after the opera- 
tion. People often balk at this, asking "How can a vector 
fixed in the structure not change when the structure is 
rotated?" The answer will be seen below. 

As a first illustration, inversion about the coordinate 
origin O in Fig. 1 simply interchanges Si and S2 (be- 
cause the spins are axial vectors; but the conclusions are 
unchanged if they are replaced by vectors). Assuming 
inversion is a symmetry of the structure one concludes 
D-SixS2 = D-S2xSi = -D-SixS2 for arbitrary Si, S2. 
Moriya's Rule 1 follows: Given this inversion symmetry, 
D = 0. (Note D does not change in these equations.) 
Next consider Rule 2. Suppose a mirror plane perpen- 
dicular to AB passes through O. Then the transformed 
spins are 

S[ = ^S2x — yS2y — ^S2z 

S'2 = ^Six — ySiy — zSiz 

Again, equating D • Si x S2 = D • S^ x S2 gives Dx = 0, 
Rule 2. This procedure can be seen to yield all 5 rules. 



under r ^ — r, making cr = d invariant under this trans- 
formation. One can see explicitly that this behavior does 
occur in the model calculation of [TO] • To accentuate the 
analogy with the DM case, we write ^ as 

p = dx(SixS2). (2) 

We now follow the same reasoning as with the energy, but 
replace it with p. That is, we consider again rotations 
that leave the structure, sites A and B plus environment, 
unchanged, and demand that p satisfy its vector prop- 
erty. We just check this against Moriya's list of (five) 
rotations, which is all possibilities that take the sites A 
and B into themselves. The result, which we found sur- 
prising, is that symmetry imposes no restrictions on a 
dipole of form ([2]). We therefore mention just one exam- 
ple to illustrate the procedure, namely inversion through 
O: The spin cross-product changes sign, d doesn't, so 
p ^ — p, which is what a vector should do; hence the 
existence of this symmetry imposes no restrictions on p. 
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FIG. 3: Spins added to Fig. 1 



FIG. 2: Highest symmetry structure with 3-fold axis along 
AB and D / 0. (D ||AB.) 

Returning to the question, "How can a vector stuck in 
the structure not rotate with the structure?" , the answer 
is that when D 7^ 0, the symmetry is so low that no 
allowed rotation will change D. As an example, we have 
constructed a structure, shown in Fig. 2, with D 7^ 0, of 
highest symmetry consistent with a 3-fold axis along AB 
(in which D || AB, according to Moriya's rule 5). (Highest 
symmetry was achieved by considering all 5 rules.) One 
sees that the full group is the 3-fold axis; there are no 
other symmetry operations. Thus D is invariant under 
the full group, i.e., it doesn^t change under any allowed 
symmetry operation. 

Let us turn from the scalar energy back to the electric 
dipole moment, a vector. Referring to ([1]), note that r 
was defined as displacement of one site relative to the 
other. For given spins Si, S2, the result p cannot depend 
on the choice of sign of r. Therefore c must change sign 



We write the other possibility for a component of the 
electric dipole moment as 

p = aSi X S2, (3) 

where a is a numerical coefficient, again a property only 
of the structure. We will go down the list of rules here, 
because the symmetry argument does give restrictions, 
depending on which operation is considered. 

1. Inversion through O 

p' = aS2 X Si = -p (4) 

again, as for the case of p under inversion, this is the cor- 
rect behavior of a vector. Therefore inversion symmetry 
imposes no restriction on p. 

2. Mirror ± AB 

The reflected spins are 

S[ = ^S2x — yS2y — ^S2z 

S2 = ^Six — ySly — zSiz 
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Using these we find 



(5) 



precisely what a vector will do under this mirror. Thus 
this ± mirror imposes no restriction on p. 
3. Mirror includes AB 

We can take this mirror as the xy plane. Then we must 
have 



^Px + yPy - ZPz 



(6) 



Since there is no interchange of Si and S2, Si x S2 trans- 
forms as an axial vector, so that 

S; X S'2 = -x(Si X S2). -y(Si X Sa)^ + z(Si x S2). (7) 

Thus this mirror symmetry implies 

p = 0. (8) 

4. 2-fold rotation axis ± AB 



Si' 



S2 



FIG. 4: Rotation of the physical elements of Fig. 3 by tt 
about the z-axis 

This axis can be taken as the z-axis, and the result of 
the rotation of the physical aspects of Fig. 3 is shown in 
Fig. 4 (for the spins projected on to the x-y plane). We 
clearly have 

s; X s'2 = -Si X S2. 

But the vector property of p demands 



d X (Si X S2), whereas it does impose restrictions on 
p = a (Si X S2). This was accomplished by following ar- 
guments closely analogous to those that lead to symmetry 
restrictions on the DM vector D. We note that just two 
of the possible symmetry operations impose a restriction 
on p, namely the mirror including AB (#3) and the 2- 
fold axis ± AB (#4). As far as we are aware, these are 
new results. Although the statement [14] that the ab- 
sence of a mirror plane perpendicular to the c-axis is an 
essential ingredient of magnetoelectricity (in connection 
with the particular material being discussed there) is a 
special case of #3, a mirror including AB. Our results 
were obtained by considering only a limited set of local 
symmetry operations, namely all those that preserve the 
pair of magnetic sites plus its environment. 

This is simpler than previous symmetry considerations 
that study restrictions on possible magnetism-caused 
electrical ordering, since they have considered the crystal 
as a whole, and consequently had to deal with the intri- 
cacies of the full crystal symmetry group and its repre- 
sentations. [H, llH, [1^ Another important advantage 
of the present approach is that it holds in any phase, 
e.g. even if the phase has been reached through a first- 
order phase transition, a situation known to invalidate 
the Landau approach. See e.g. fi^ . 

While our arguments about p and p are closely analo- 
gous to those leading to information about D, one should 
recognize an important difference, namely the electric 
dipole is a local physical variable, the DM D being a 
parameter in the spin Hamiltonian. Thus, considering 
a given crystal structure, a question with no meaning 
concerning D, but with serious meaning concerning the 
dipole, is what kind of ordering, or spontaneous symme- 
try breaking can occur among the macroscopic number of 
dipoles in the crystal. An important illustration is to con- 
sider translating Fig. 1 to form a linear chain, with the 
corresponding extension of Fig. 3, but rotating the spins 
to form a long-range ordered spiral. The contribution to 
the dipole moment from every pair is identical, for either 
p or p, so that both could contribute to a macroscopic 
electric polarization P [gl, 0, [ij]. This, extended to 3 
dimensions, is one of the bases for understanding much 
of the recently observed multiferroic behavior connected 
to spiral-like magnetism. 



p = -Xp^ -yPy +Zp^. 

Conclusion: Pz = 0. That is, this 2-fold rotation symme- 
try axis demands p± the 2-fold axis. 
5. AB is an n-fold axis, n > 2 

This imposes no restriction on p: There is no interchange 
of Si and S2. Thus Si x S2 transforms as a pseudovec- 
tor which, under this proper rotation, is the same as a 
vector. 

These results were checked against the simplest of the 
microscopic models in [10], namely 2 hydrogen atoms, 
convenient because of the simplicity. 

To summarize, we have shown quite rigorously and 
generally that symmetry imposes no restrictions on p = 



^ 1 X 



FIG. 5: Triangular plaquette of spins (darker arrows) and 
electric dipoles (lighter arrows). 

We conclude by mentioning some other physical con- 
sequences of our basic results. 
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(a) The robustness of p under symmetry requirements 
may be why it has been found experimentahy in many 
different materials, whereas p has been found, as far 
as we are aware, only in one, namely RbFe(MoO)4)2 
(RFMO) [14]. 

(b) In connection with RFMO (also discussed in [15]): 
At low temperature it contains triangular lattice planes 
of spins, the spins lie in these planes and form a 120°- 
turn- angle spiral ("the 120° structure"). For the crystal 
structure see [20|; the low-temperature space group is P3. 
Further, there is no mirror plane parallel to these planes, 
and no 2-fold axis, so according to our result. Rules 3 
and 4, there can be a local electric dipole moment p, 
which lies ± these planes (since the spins lie in these 
planes and by definition of p). Also, since the handed- 
ness of the cross-product is coherent, i.e. each nearest 
neighbor spin pair gives the same dipole moment, the 
observed ferroelectric macroscopic polarization, P ± to 
the planes, is implied. The authors note [M] that the 
existence of a 3-fold axis ± to the planes (the c-axis) im- 
plies there can not be a component of P parallel to the 
planes. On the other hand, the robustness of p suggests 
strongly that there are local electric dipoles that order 
antiferroelectrically in the 120° state as indicated for a 
single triangular plaquette in Fig. 5. The location of the 
electric moments at the midpoints of the triangle edges 
is symbolic of the actual bond charge density found in 
the microscopic models of 

U K im (although, with the 
exception of [10] b, the high symmetry assumed in these 
calculations requires p = 0). Such a charge distribution 
would be ordered in the crystal (it's tied strongly to the 
magnetism), and therefore detectable, perhaps even in 
practice. 

These last remarks are predicated on the view that if a 
physical effect is allowed by symmetry, then it will occur 



in nature. This connects to the question of the micro- 
scopic origin of the observed effects. The microscopic 
theories of [H, [13, [lH fundamentally consider a single 
pair of magnetic atoms; and the mechanism for obtain- 
ing an electric dipole moment is the quantum-mechanical 
calculation of the charge density treating the atomic elec- 
trons as localized, and considering the spin-orbit coupling 
as a perturbation. As we noted, an essential difference 
between the work of 0, [ij] and that of 0] is in the 
role of the spin-orbit coupling: in the former case, this 
contributes to matrix elements between Wannier func- 
tions localized (mainly) on one atom (we have called this 
"intra-atomic SO coupling"), whereas in the latter, the 
SO coupling is the operator in the interatomic hopping 
integral ("inter-atomic SO coupling"). One can see anal- 
ogous contributions to the DM D il8j. As we have noted, 
in the earlier calculations [H, |Tl|, ^^a, the over- all sym- 
metry was assumed to be that of two dots in free space, 
and this led to A = p; when lower symmetry was consid- 
ered later in the inter- atomic SO case [10] b, the contri- 
bution p to A was found to occur. Due to the generality 
and rigor of the new symmetry theory presented here, 
there is little doubt that consideration of lower symme- 
try in the intra-atomic SO coupling case will also lead to 
a contribution of the form p. We note that the symmetry 
effects show up in our calculation [10] via the potential 
energy V^(r) that appears in the SO operator W x p • s. 
We have added these comments because the intra-atomic 
model appears to give a quantitatively much larger elec- 
tric dipole moment [11], and we would hope that that 
approach will be pursued considering lower symmetry. 

We thank Mr. Z. Rak for help in understanding the 
RFMO structure. Helpful communications with A. B. 
Harris, G. Lawes, M. Kenzelmann, Y. Tokura, N. Na- 
gaosa, and H. Katsura are gratefully acknowledged. 
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